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Abstract. Using a change-of- measure argument, we prove an equality in law be- 
tween the process of largest eigenvalues in a generalized Wishart random-matrix 
process and a last-passage percolation process. This equality in law was conjectured 
by Borodin and Peche (2008). 



1. Introduction 

The past decade has witnessed a surge of interest in connections between random 
matrices on the one hand and applications to growth models, queueing systems, 
and last-passage percolation models on the other hand; standard references are 
Baryshnikov (2001) and Johansson (2000). In this note we prove a result of this 
kind: an equality in law between a process of largest eigenvalues for a family of 
Wishart random matrices and a process of directed last-passage percolation times. 

To formulate the main result, we construct two infinite arrays of random vari- 
ables on an underlying measurable space, along with a family {P 7r,7r } of probabil- 
ity measures parametrized by a positive ./V-vector tt and a nonnegative sequence 
{fr n : n > 1}. The elements of the first array {Ay : 1 < i < N, j > 1} are in- 
dependent and Aij has a complex zero-mean Gaussian distribution with variance 
l/(iTi + 7Tj) under P 7r > 7r . That is, both the real and complex part of Ay have zero 
mean and variance l/(2ni+2'Kj). Write A(n) for the Nxn matrix formed by the first 
n columns of A, and define the matrix-valued stochastic process {M(n) : n > 0} 
by setting M(n) = A(n)A(n)* for n > 1 and by letting M(0) be the N x N zero 
matrix. We call {M(n) : n > 0} a generalized Wishart random-matrix process, 
since the marginals have a Wishart distribution if tt and tt are identically one and 
zero, respectively. 

The elements of the second array {Wij : 1 < i < N,j > 1} are independent and 
Wij is exponentially distributed with parameter 71".; + itj under p 7r ' 7r . We define 

Y(N, n) = max Wa, 
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where Yl(N,n) is the set of up-right paths from (1,1) to (N,n). The quantity 
Y(N, n) arises in last-passage percolation models as well as in series Jackson net- 
works in queueing theory, see for instance Dieker and Warren (2008) or Johansson 
(2009). 

The following theorem, a process-level equality in law between the largest eigen- 
value of M(n) and Y(N,n), is the main result of this note. Given a matrix C, we 
write sp(C) for its vector of eigenvalues, ordered decreasingly. 

Theorem 1.1. For any strictly positive vector ir and any nonnegative sequence w, 
the processes {sp(M(n))i : n > 1} and {Y(N,n) : n > 1} have the same distribution 
under P 7r > 7r . 

It is known from Defosscux (2008); Forrester and Rains (2006) that this holds 
in the 'standard' case, i.e., under the measure P :— pC 1 '---' 1 )^ ' '---). i n its stated 
generality the theorem was conjectured by Borodin and Peche (2008), who prove 
that the laws of Y(N, n) and the largest eigenvalue of M(n) coincide for fixed n > 1. 
Our proof is based on a change-of- measure argument, which is potentially useful to 
prove related equalities in law. 

Throughout, we use the following notation. Wc let Hn,n be the space of all 
N x N Hermitian matrices, and W N the set {x £ M. N : x\ > ... > x^}. For 
x, x' £ W N , we write x -< x' to mean that x and x' interlace in the sense that 



2. Preliminaries 

This section provides some background on generalized Wishart random matri- 
ces, and introduces a Markov chain which plays an important role in the proof of 
Theorem 1.1. 

2.1. The generalized Wishart random-matrix process. Under P 1 ' 1 , the generalized 
Wishart process {M(n) : n > 0} from the introduction has independent increments 
since, for m > 1, 



where Aj m is the complex conjugate of Aj m . In particular, the matrix- valued 
increment has unit rank. The matrix M(m) — M(m— 1) can be parameterized by 
its diagonal elements together with the complex arguments of Ai, n for 1 < i < N; 
under P^^ , these are independent and the former have exponential distributions 
while the latter have uniform distributions on [0, 2-k\. (This fact is widely used in 
the Box-Muller method for computer generation of random variables with a normal 
distribution.) Since the i-th diagonal element has an exponential distribution under 
P^^ with parameter i^i + 7r m , we obtain the following proposition. 

Proposition 2.1. For any m > 1, the P v ' w -law of M(m) — M(m— 1) is absolutely 
continuous with respect to the P-law of M(m) — M(m — 1), and the Radon-Nikodym 
derivative is 



x[ > X\ > x' 2 > X2 > ■ • • > x' N > Xn- 



M(m) = M(m - 1) + (A im A f 



l 3 m )l<ij<N ' 



(2.1) 
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2.2. A Markov transition kernel. We next introduce a time-inhomogeneous Markov 
transition kernel on W N . We shall prove in Section 3 that this kernel describes 
the eigenvalue-process of the generalized Wishart random-matrix process of the 
previous subsection. 

In the standard case (ir = 1, n = 0), it follows from unitary invariance (see 
Dcfosscux (2008, Sec. 5) or Forrester and Rains (2006)) that the process {sp(M(n) : 
n > 0} is a homogeneous Markov chain. Its one-step transition kernel Q(z, •) is 
the law of sp(diag(z) + G), where G = {gi(jj : 1 < i,j < N} is a rank one matrix 
determined by an iV-vector g of standard complex Gaussian random variables. For 
z in the interior of W N , Q(z,-) is absolutely continuous with respect to Lebesgue 
measure on W N and can be written explicitly as in Defosseux (2008, Prop. 4.8): 

Q^dz') = ^le-^'^h^ydz', 

where A(z) :— Yli<i<j<N( z i — z j) i s tne Vandermonde determinant. 

We use the Markov kernel Q to define the aforementioned time-inhomogeneous 
Markov kernels, which arise from the generalized Wishart random-matrix process. 
For general tt and n, we define the inhomogeneous transition probabilities Q^-i „ 
via 

Ql^Jz, dz>) = Ufa + n n )^le-^ X£M-»)Q {Z , dz'), 

dct{ e -^Q 

K{Z) - A(n)A(z) ■ (2 - 2) 
Note that h^{z) extends to a continuous function on (0, oo)^ x W N (this can im- 
mediately be seen as a consequence of the Harish-Chandra-Itzykson-Zuber formula, 
see (3.2) below). 

One can verify that the Qf^ are true Markov kernels by writing 1{ Z ^ Z /} = 
det{l{ 2i<z ^}} and applying the Cauchy-Binet formula 

det{£i(zj)} det{tp j(zi)}dz = det I / ^(z)ip j (z)dzl . 
iw N Urn. J 

3. The generalized Wishart eigenvalue-process 

In this section, we determine the law of the eigenvalue-process of generalized 
Wishart random-matrix process. Although it is not essential to the proof of The- 
orem 1.1, we formulate our results in a setting where sp(M(0)) is allowed to be 
nonzero. 

Write to m for the 'uniform distribution' on the set {M E Hn.n ■ sp(M) = 
fi}. That is, to m is the unique probability measure invariant under conjugation by 
unitary matrices, or equivalcntly is the law of U diag(/i)[/* where U is unitary 
and distributed according to (normalized) Haar measure. We define measures 
by letting the P^<*-law of {M(n) - M(0) : n > 0} be equal to the P^-law of 
{M(n) : n > 0}, and letting the P^-distribution of M(0) be independent of 
{M(n) — M(0) : n > 0} and absolutely continuous with respect to to m with Radon- 
Nikodym derivative 

_£^_ e - Ef =1 * eX p(-tr[(diag(7r) - 7)M(0)]), (3.1) 
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where cat is a constant depending only on the dimension TV and / is the identity 
matrix. Recall that h^(p) is defined in (2.2). That this defines the density of a 
probability measure for all it and /x follows immediately from the Harish-Chandra- 
Itzykson-Zuber formula (e.g., Mehta (2004, App. A. 5)) 

/ exp (-tr(diag(7r)J7diag(/x)?7*)) dU = c^XC"), (3.2) 
Ju 

writing dU for normalized Haar measure on the unitary group. Throughout, we 
abbreviate pU.-. 1 )^ . .-) by p ^ Notc that thc P ^*_\ aw an d the P^-law of 
{M(n) : n > 0} coincide if fi = 0. 

The following theorem specifies the P^'Maw of {sp(M(n)) : n > 0}. 

Theorem 3.1. For any fi G W N , {sp(M(n)) : n > 0} is an inhomogeneous 
Markov chain on W N under P^ ,7r ' , and it has the Qn-i, n °f Section 2.2 for its 
one-step transition kernels. 

Proof. Fix some /x e W N . The key ingredient in the proof is a change of measure 
argument. We know from Dcfosseux (2008) or Forrester and Rains (2006) that 
Theorem 3.1 holds for the 'standard' case 7r = (1, . . . , 1), it = 0. 

Writing P*>* and P n for thc distribution of (M(0), . . . ,M(n)) under P^<* and 
P M respectively, we obtain from Section 2.1 that for n > 0, 

^(M(0),...,M(n)) 
MaO 

x exp ^-tr((diag(7r) - I)M(n)) - ^ it m tr(M(m) ~ M(m - 1))^ , 

where C Vt %(n,N) = YliLi Ilj=i( 7r i + ^j)- Let thc measure p 7 ^ 71 (and p n ) be the 
restriction of P^ % (and P„) to the cr-field generated by (sp(M(0)), . . . , sp(M (ra))). 
Then we obtain for n > 0, 

-(sp(M(0)),...,sp(M(n))) 



E 



Pu 



dPr, 



(M(0),...,M(n)) 



s P (M(0)),...,sp(M(n)) 



where Ep^ denotes the expectation operator with respect to P^. Since the P M - 
distribution of (M(0), . . . , M{n)) given the spectra is invariant under component- 
wise conjugation by a unitary matrix U, we have for /j = ~< jjS 1 ^ -<...-< 



E F 



exp (-tr(diag(Tr)M (n))) sp(M (0)) = A« (0) , • • ■ , sp(M (n)) = xt (n) 
exp (-tr(diag(7r)J7diag(At (n) )?7*)) dU 
= c^M/i (n) ), 

where the second equality is the Harish-Chandra-Rzykson-Zuber formula. From 
the preceding three displays in conjunction with tr(M) = sp(M)j, we conclude 
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that 

*-Vm (1) ,...,m (b) ) 



dp n 

N n N 



= C7r ' #(n,7V) exp ( _ £5> N r) -^"1 +5> (B) ■ 

Since sp(M(-)) is a Markov chain with transition kernel Q under P^, we have 
P;- # (sp(M(l)) G d^ l \ . . . , sp(M(n)) G d^ (n) ) 

= . . . , M H)P AJ ( sp (M(l)) G d// 1 ), . . • ,sp(M(n)) G 

the last equality being a consequence of the definition of Q^-x k and the expression 
for dp^ % /dp n . □ 

4. Robinson-Schensted-Knuth and the proof of Theorem 1.1 

This section explains the connection between the infinite array {Wy} of the 
introduction and the Markov kernels Qn-i, n - ^ n conjunction with Theorem 3.1, 
these connections allow us to prove Theorem 1.1. 

The RSK algorithm. The results in this section rely on a combinatorial mechanism 
known as the Robinson-Schensted-Knuth (RSK) algorithm. This algorithm gener- 
ates from a p x q matrix with nonnegative entries a triangular array x = \x\ : 1 < 
j < Pi 1 < i < j} called a Gelfand-Tsetlin (GT) pattern. A GT pattern with p 
levels x 1 , . . . , x p is an array for which the coordinates satisfy the inequalities 

x k x k-l — X k-1 — x k-2 — ■ ■ ■ — x 2 ^ X l X l 

for fc = 2,...,p. If the elements of the matrix are integers, then a GT pattern 
can be identified with a so-called semistandard Young tableau, and the bottom row 
x p = {x?; 1 < i < p} of the GT pattern corresponds to the shape of the Young 
tableau. We write K p for the space of all GT patterns x with p levels. 

By applying the RSK algorithm with row insertion to an infinite array : 
1 < i < N, 1 < j < n} for n = 1,2, . . ., we obtain a sequence of GT patterns 
x(l), x(2), .... It follows from properties of RSK that 

x i(n) = max (4.1) 
(ij)ep 

where H(N,n) is the set of up-right paths from (1, 1) to (N,n) as before. Details 
can be found in, e.g., Johansson (2000) or Dieker and Warren (2008, case A). 

Greene's theorem generalizes (4.1), and gives similar expressions for each compo- 
nent of the pattern x' i (n), see for instance Chapter 3 of Fulton (1997) or Equation 
(16) in Doumerc (2003). As a consequence of these, we can consider the RSK algo- 
rithm for real-valued &j and each x(n) is then a continuous function of the input 
data. 
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We remark that the RSK algorithm can also be started from a given initial GT 
pattern x(0). If RSK is started from the null pattern, it reduces to the standard 
algorithm and we set x(0) = 0. 

The bijective property of RSK. RSK has a bijective property which has important 
probabilistic consequences for the sequence of GT patterns constructed from spe- 
cially chosen random infinite arrays. Indeed, suppose that {£y : 1 < i < N, j > 1} 
is a family of independent random variables with £y having a geometric distribu- 
tion on Z + with parameter aibj, where {ai : 1 < i < N} and {bj : j > 1} are two 
sequences taking values in (0,1]. Write {X(n) : n > 0} for the sequence of GT 
patterns constructed from £. 

Using the bijective property of RSK it can be verified that the bottom rows 
{X N (n) : n > 0} of the GT patterns evolve as an inhomogeneous Markov chain 
with transition probabilities 

i=1 s x (a) 
where s\(a) is the Schur polynomial corresponding to a partition A: 

s x( a ) = flX ' 
with the weight a x of a GT pattern x being defined as 

N 

«-=^n-F' f " E " f_1 - 

fe=2 

This is proved in O'Connell (2003) in the special case with bj = 1 for all j, and the 
argument extends straightforwardly; see also Forrester and Nagao (2008). 

Non-null initial GT patterns generally do not give rise to Markovian bottom-row 
processes. Still, the inhomogeneous Markov chain of bottom rows can be con- 
structed starting from a given initial partition A with at most N parts by choosing 
X(0) suitably from the space of a GT patterns with bottom row A: X(0) should be 
independent of the family } with probability mass function 

a x 

s\(a) 

Exponentially distributed input data. We now consider the sequence of GT patterns 
{X^(n) : n > 0} arising from setting <Zj = 1 — TTi/L and bj = 1 — nj/L in the above 
setup, and we study the regime L — > oo after rescaling suitably. In the regime L — > 
oo, the input variables {£,ij/L} (jointly) converge in distribution to independent 
exponential random variables, the variable corresponding to £y /L having parameter 
7Tj +nj. Thus, the law of the input array £ converges weakly to the P^—law of the 
array {Wij ■ 1 < i < N,j > 1} from the introduction. Refer to Doumerc (2003) 
and Johansson (2000) for related results on this regime. 

By the aforementioned continuity of the RSK algorithm and the continuous- 
mapping theorem, {Xt(n)/L : n > 0} converges in distribution to a process {Z(n) : 
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n > 0} taking values in GT patterns with TV levels. As a consequence of the above 
results in a discrete-space setting, we get from (4.1) that 

Z?(n) = max V W M . 

(lj)£P 

Moreover, the process of bottom rows {Z N (n) : n > 0} is an inhomogeneous Markov 
chain for which its transition mechanism can be found by letting L — > oo in (4.2): 

Lemma 4.1. Under P Tr > 7r 7 the process {Z N (n) : n > 0} is an inhomogeneous 
Markov chain on W N , and it has the Qn-i n of Section 2.2 for its one-step transi- 
tion kernels. 

A similar result can be obtained given a non-null initial bottom row /i G W N . 
In case the components of \x are distinct, the distribution of the initial pattern 
Z(0) should then be absolutely continuous with respect to Lebesgue measure on 
{z g Kat : z N — fi} with density 

detje-^w} ' 

where c = (e _7ri , . . . ,e~ 7TN ). 

Proof of Theorem 1.1. We now have all ingredients to prove Theorem 1.1. We 
already noted that Z^ (n) equals Y(N,n) . Thus, for any strictly positive vector it 
and any nonnegative sequence it, {Y(N,n) : n > 1} has the same P^'^-distribution 
as {Zi(n) : n > 1}. In view of Theorem 3.1 and Lemma 4.1, in turn this has the 
same P^-distribution as the largest-eigenvalue process {sp(M(n))i : n > 1}. This 
proves Theorem 1.1. 



References 



Yu. Baryshnikov. GUEs and queues. Probab. Theory Related Fields 119, 256-274 (2001). 
A. Borodin and S. Peche. Airy kernel with two sets of parameters in directed percolation 

and random matrix theory. J. Stat. Phys. 132, 275-290 (2008). 
M. Defosscux. Orbit measures, random matrix theory and interlaced determinantal pro- 
cesses (2008). ArXiv:0810.1011. 
A. B. Dieker and J. Warren. Determinantal transition kernels for some interacting particles 

on the line. Ann. Inst. H. Poincare Probab. Statist. 44, 1162-1172 (2008). 
Y. Doumerc. A note on representations of eigenvalues of classical Gaussian matrices. 

In Seminaire de Probabilites XXXVII, volume 1832 of Lecture Notes in Math., pages 

370-384. Springer, Berlin (2003). 
P. J. Forrester and T. Nagao. Determinantal correlations for classical projection processes 

(2008). ArXiv:0801.0100. 
P. J. Forrester and E. M. Rains. Jacobians and rank 1 perturbations relating to unitary 

Hessenberg matrices. Int. Math. Res. Not. page Art. ID 48306 (2006). 
W. Fulton. Young tableaux. Cambridge University Press, Cambridge (1997). 
K. Johansson. Shape fluctuations and random matrices. Comm. Math. Phys. 209, 437-476 

(2000). 

K. Johansson. A multi-dimensional Markov chain and the Meixner ensemble. Arkiv for 

Matematik, to appear (2009). Doi:10.1007/sll512-008-0089-6. 
M. L. Mehta. Random matrices. Elsevier/ Academic Press, Amsterdam, third edition 

(2004). ISBN 0-12-088409-7. 
N. O'Connell. Conditioned random walks and the RSK correspondence. J. Phys. A 36, 

3049-3066 (2003). 



